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Abstract. 

We present a method that yields three decoupled covariant equations for three 
complex scalars, which completely govern electromagnetic perturbations of non- 
vacuum, locally rotationally symmetric class II spacetimes. One of these equations 
is equivalent to the previously established generalized Regge- Wheeler equation for 
electromagnetic fields. The remaining two equations are a direct generalization of the 
Bardeen-Press equations. The approach undertaken makes use of the well established 
3+1 (and 2+1+1) formalism, and therefore, it is an ideal setting for specifying 
interpretable energy-momentum on an initial spacelike three-slice as the perturbation 
sources to the resultant electromagnetic radiation. 
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1. Introduction 

In this paper, we consider electromagnetic (EM) perturbations to non-vacuum 
background spacetimes that exhibit local rotational symmetry of class II [ol IT2*1 IT4*1 I2U|. 

The 3+1 slicing formalism^] is employed as it has the desirable property of 
retaining a transparent and physical interpretation of the energy-momentum tensor. 
This corresponds to the "zero rotation case" for the 1+3 threading formalism expounded 
by JB] • The four- dimensional spacetime is foliated into a family of spacelike three-slices, 
and subsequently, all tensors are decomposed irreducibly into spatial parts residing 
within the three-slices, and parts which are orthogonal. As the assumption that 
the background is locally rotationally symmetric (LRS) class II is imposed, there is 
additional structure hidden within the equations, which is difficult to exploit in 3+1 
form. To resolve this, a further splitting of the 3+1 quantities is desirable. Every 
three-slice is therefore foliated by a family of two-slices, and all spatial quantities are 
decomposed into "angular" parts residing entirely on the two-slices and "radial" parts 
orthogonal to the two-slices. Analogous to the difference between the 3+1 and 1+3 
formalisms stated above, this will be referred to as a 2+1+1 slicing formalism, and this 
is a special case of the 1+1+2 threading formalism developed in [HI CHj. Ultimately, 
the LRS class II equations reduce to a system of partial differential equations (PDEs) 
involving only covariant scalar quantities^- 

First we use linear algebra techniques to show that Maxwell's equations decouple 
naturally by constructing complex quantities. The 2+1+1 splitting is then used to 
decoupling the "radial" quantities from the "angular" quantities. The "radial" quantity 
satisfies the Regge- Wheeler (RW) equation [TH] generalized to LRS class II spacetimes, 
for which the real and imaginary equations correspond to those derived in [Bj. 

This paper extends the work in [B] by performing an additional covariant splitting 
of the all "angular" quantities residing on the two-slices. Two covariant, decoupled, 
angular equations are ultimately derived, and we shown these are a direct generalization 
of the Bardeen-Press (BP) equations [2j toward non-vacuum LRS class II spacetimes. 

1.1. Mathematical Preliminaries 

Greek indices run from to 3, lower case Latin from 1 to 3 and capital Latin from 2 to 3. 
The standard symmetric ( ) and skew-symmetric [ ] notation is used, and angle brackets 
<> represent the spatially projected onto a three-surface, symmetric, and trace- free part 
of a tensor. We define the four-coordinates x^ := (t,r, x A ), where t is the "temporal" 
coordinate, r is the "radial" coordinate and x A are the "angular" coordinates. Note 
that, temporal, radial and angular are used generically here. Geometric units are used 
whereby 8ttG = c = 1. The completely anti-symmetric Levi-Civita pseudo-tensor (e Mi/crr ) 
is defined such that e a/ 3 7W e a/37a; = —4!. The Lie derivative along a vector field 
is denoted Cx- The complex number is denoted by i. The covariant derivatives with 
respect to the full spacetime, three-surface and two-surface are respectively (V M , D i: c^)- 
The Riemann tensors are therefore 2V[ M V„]V^ = R f _ iUaT V T , 2 D^D^Wa = (3> R i _ LV(T tW t 
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and 2d^d v \V a = R^ V otV t where W a has been projected onto the three-surface, and 
V a has been projected onto the two-surface. The Bianchi identities for a two-surface 
imply that the two-Ricci tensor can be written in terms of the Gaussian curvature (K g ) 
according to jR.^ = K g S^ u where represents the metric of the two-surface. 

2. The Background LRS Class II Spacetime 

The background spacetime is assumed to be LRS class II. This section first presents 
the non-trivial equations for the 3+1 formalism and subsequently, the 2+1+1 system of 
equations is discussed. 



2.1. 3+1 Formalism 

The comprehensive 3+1 formalism for general relativity consists of both the Einstein 
field equations and the Bianchi identities. The 3+1 formalism is well established 
throughout the literature; for example, see [U El ESj and the "zero rotation cases" of 
the 1+3 formalism in Thus only the essential mathematical tools required 

to decompose equations into 3+1 form are presented here, followed by a summary of 
the equations. The construction of spatial three-surfaces requires a time-like normal 
defined, 

= —aV^t and n a n a := —1, (1) 

where a > is the lapse function and the Cauchy time function has, without loss of 
generality, been set equal to t. The projection tensor is used to project any quantity 
onto the three-surfaces and is defined, 

J- Mi ,:= g^u + n^n u . (2) 

The associated three-metric permitted for a spatial three-slice is Ly. The extrinsic 
curvature of the three-surfaces is related to the covariant derivative of the normal, which 
is decomposed according to, 

= —A^ v - i K - n^h v . (3) 

Here, A^ u and K are the trace-free and trace parts of the extrinsic curvature, and 
ri M := n a V (1 n (J is the four-acceleration. The ten independent components of the energy- 
momentum tensor (T M!/ ) are defined 



p := T Q/3 r/V, P := i ± a(3 T a/3 , 

j[i '■= —-L f i /3 Tp a n a , 7T^ U := + M a + l/ /3 T Q , /3 — - + Mi ,+ a/3 T a p. 



(4) 



3 

Here p is the mass-energy density, P is the isotropic pressure, is the spatial mass- 
energy flux and TT^ represents the spatial, trace-free, anisotropic pressure. 

The ADM equations P comprise of both conservation of energy-momentum, and 
the Einstein field equations, 

V% a = and R fll/ --g^R + Ag fll/ = T flv , (5) 
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where A is the cosmological constant. These can be decomposed into their 3+1 form 
and are respectively, the continuity and Euler equations, 

(£ n - K) p + (D k + 2 h k ) Jk - 7T km A km - P K = 0, (6) 

(C n - K)ji + {D k + h k )7T ik + (A + h i )P + ph l = 0. (7) 

The Einstein field equations consist of the Hamiltonian and momentum constraints, and 
evolution equations for the extrinsic curvature, 

1 (3) 1 1 
P= 2 R ~ 2 AkmAkm + 3 R2 ~ A > ^ 

3l = D k (A lk -^UkK), (9) 
P = - l -K)K- l -A km A km + ^ (D k + h k )h k -\R + K (10) 

1 (3) 

TTij = -{£>n - 2 K ) A <V> - ( D <i + n <i) n 'j> + Qij- ( U ) 

Here, ^Qij is the trace-free part of the three- Ricci tensor on a three-slice and <3 ^2 is the 
corresponding three-Ricci scalar. The evolution equation for the three-metric is given 
by 

2 

l~"rt -Ly ^ Aij — _Ljj K , (12) 

and additional equations of state for the energy-momentum are essential to close this 
system in general. 

The ADM equations, given by (P)- (fT2"|) . can be supplemented with the 3+1 gravito- 
electromagnetic (GEM) formalism [HO HE]. The Weyl tensor (C^ uaT ) is split into two 
GEM tensors defined according to 

Eftv := C a ^ u n a n p and B^ v := -e^C^ uv n u , (13) 

where the completely anti-symmetric three-Levi-Civita tensor can be been defined to 
satisfy e^ ua := e a ^ U(T n a . The assumption of LRS class II spacetimes is sufficient to 
enforce the gravito-magnetic tensor (B^) to vanish everywhere [14j. Therefore, using 
LRS class II symmetry, it follows that the non-trivial 3+1 GEM system^ reduces to 

D k [E ki + ±7T k ^ = - l -e kim A k n (e nmp Jp ) + ±D iP + ^ K, (14) 
[c n - ^K^j E Kij> + ^e km<l (D k + 2n k ) ( Cj 

+5A <t k ^E j>k + ^7T j>k J + \ (c n + ^Kj Tl <l3> = \(P + p)A r (15) 

Finally, there are some additional relationships which yield useful information regarding 
the coupling between the ADM and GEM systems§, 

Eij — —71 ' ij = CnA^jy + A k< iAjy + (-D<j + h < i)hj > , (16) 

| The 3+1 GEM system is the (once contracted) Bianchi identities decomposed into 3+1 form. 

§ 1)16(1 and 1)17(1 are found by first expressing the definition for the gravito-electric tensor (given in 

in terms of the Riemann tensor, Ricci tensor and Ricci scalar. Subsequently, this is decomposed to 



mn-\ 
C J> Jn) 
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^ij + 2^"^ = ~ ^-k<i^j> k + T^AijK. (17) 

The assumption of an LRS class II spacetime will greatly simplify these equations as 
most of them will have many trivial components, and this can certainly be analyzed from 
a component point of view. However, the covariant 2+1+1 formalism will be responsible 
for revealing only non-trivial scalar equations and this leads to simplifications. 



2.2. 2+1+1 Decomposition 

In LRS class II spacetimes, every three-slice can be further foliated by a family of two- 
slices. A normal to these two-slices (iV 7 *) is spatial (i.e. N a n a = 0), and is normalized 
according to, 

± af3 N a Np : = 1, (18) 
from which it follows immediately, 

N»=-±=D»r, (19) 



where V + 11 > and r is the "radial" coordinate]). A projection tensor which projects 
quantities onto the two-surfaces is therefore 

S^ V :=^ U -N^N V . (20) 

The extrinsic curvature of the two-surfaces is related to the covariant three-derivative 
of the spatial normal, 

D^N V = -A^ - i<S M „ K. + ivyi, (21) 

where A^ and K are the trace-free and trace parts of the extrinsic curvature of the 
two-surface, and := iV"^^. 

The 3+1 equations can now be decomposed into 2+1+1 form. The trace-free 
extrinsic curvature of the three-slice, the gravito-electric tensor and the anisotropic 
stress are decomposed as, 

= ^~\ + 2 Y, {lx N v) + S N„N V , (22) 
= Ep,-\ S ^ E + 2 e b N ») + 8 N » N »> ( 23 ) 

iffiu = - x - s^n + 2 + n n»n v , (24) 

where £ := A aP N a N?, £ M := S^A a pN^ = S «S v p A a p + E, and analogous 
definitions follow for E^ v and TT^. It is important to note that now E^, £^ v and IT^ 
are trace-free with respect to the two-slice; for example, S a ^H a p = 0. 

give an evolution equation for the extrinsic curvature. Finally, linear combinations with Ijlljl gives the 
required results. 

|| The general solution for this normal is = f D^gif) where / = f(t,r,x A ) and / g'(r) = 1/V-L 11 - 
This is a similar situation to the Cauchy time function; g(r) must be a function of r only and can, 
without loss of generality, be set to equal r. 
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The general 2+1+1 equations are not presented here, as only the LRS class II 
equations are required. LRS class II implies that the only non-vanishing terms are 
scalar quantities jH], i.e. 

LRS class II: {S, K, A, 8, K, A, p, P, J, II}, (25) 

where A := h a N a and J := j a N a . A more comprehensive 1+1+2 decomposition can 
be found in jOHOj- Following the work of it is instructive to work with the Gaussian 
curvature of the two-surface, 

K g = -8 + \K} - (\k - hlf + \(p + A) - in. (26) 
Therefore, the Ricci curvature for the three-surface can be decomposed according to 

(3) R = 2 (c N - -ICj )C + 2K g , (27) 

= -\s^(C N K -2K g ) + \n^N v (C n 1C -2 K g ). (28) 
6 3 

In summary, the non-trivial LRS class II equations after irreducible decomposition are 

the continuity and Euler equations 

(C n - S - i K)J + (C N + A)P + (C N + A - | /C)n + P A = 0, (29) 
(C n - K)p + (C N + 2 A - K)J - K P -^EU = 0. (30) 



The Hamiltonian and momentum constraints, the Raychaudhuri equation [TH] and an 
evolution equation for E are respectively, 

C N -\k\k. + K g + i/T 2 - i £ 2 = p + A, (31) 

C N -^)CjZ-^C N K = J, (32) 

(Cn - l -K)K +(C N + A- K)A -h?= l -(p + W- 2A), (33) 

-(C n + ~E)(E - 2 -K) + (C N -A- ~K)K = U + p + P (34) 
The constraint and evolution equations governing the gravito-electric tensor are 

(C N - Ik)(£ + in) - i^ATP = -i(E - ~*Q^ (35) 

(A. + ^E - K)£ - i£„p + i(£„ + is - itf)II = -i/CJ + i( p + P)(S - j^). (36) 

The equations giving information regarding the coupling of the ADM and GEM systems 
reduce to 

8 - in = (C n - l -K + is)(S - 2 -K) + AK. + i(p + 3 P - 2 A), (37) 

5 + in = - i/C)/C - (iff + S)(S - |l0 - \{ P + A). (38) 
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Finally, an additional evolution equation for /C is^f, 

Un - \k + ~e\ K + (e - 2 -K\ A = -J. (39) 

These are the 2+1+1 equations governing the dynamics of the background spacetime" 1 ". 
They need to be prescribed in order to evaluate the perturbations. 

2. 3. Commutators 

The decoupling of Maxwell's equations relies heavily on the commutation relationships 
between the Lie derivatives along both and iV M , and the covariant two-derivative on 



the two-surface. For any scalar function /, 

(C N + A)C n f - (C n - E - i K)C N f = 0, (40) 

d a C n f - C n d p f = 0, (41) 

d a C N f - C N d p f = 0, (42) 

dfrd v] f = 0. (43) 

Here, the "bar" is used to denote an index which is projected onto the two-surface, i.e. 
Vp := Sfj, a V a . Similarly, for any two-tensor (i.e. N a <fr a = n a & a = 0), 

[c N + A + I/choc - h: + \k)^\ = [c n - 1 E]^ + (44) 

d a C n $ D - C n d-^u = 0, (45) 

d^C N <$> D - C N dj& p = 0. (46) 



3. Maxwell's Equations 

There are two primary ideas discussed in this section. The first discusses Maxwell's field 
equations from a general point of view and illustrates the first stage of decoupling. The 
latter then focuses on the perturbed quantities and deriving the generalized RW and 
BP equations. 

3.1. Maxwell's Equations in General 

The covariant 3+1 formalism for EM fields has been studied previously[7J 122 I2H|- 
The spatial EM field intensities (E^ and B^) are expressed covariantly in terms of the 
anti-symmetric Faraday tensor (F^), 

E„:=F m n a and B» := \e m ^ . (47) 
It then follows that the EM tensor is decomposed into 3+1 form according to 

F uu = e uua B a -2E {ll n u] . (48) 

% This may be derived by decomposing n M <S l " T (2V[ A1 V i ,]-/V CT — R^ v <7tN t ) — as indicated in 0. 

+ These equations correspond to those presented in [HI E] where the transformation from our notation 

to theirs is {£, if, K., p, J, n^, N^} — > {— £, -if, —<f>, fi, Q,u M , n^}. 
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Maxwell's equations are expressed conveniently in terms of F^ v through 

J... 

2 



V [ff F H =0 and V a F m = (49) 



Here J M represents the charge- current four-vector which is decomposed as 

J M =_L J M - {n a J a )nn = v + en M , (50) 

where _L J M := i^ is the three-current density and n a J a := —e is the charge density. 
Maxwell's equations can now be expressed in 3+1 form, 

D k B k = 0, (51) 

D k E k = i e , (52) 

(C n -lK)B l + e^(D k + n k) E m + 2A*B k = 0, (53) 

(£ n - ijC)^ - e km (D k + h k )B m + 2 A k E k = ~ U. (54) 

This is a coupled system of two constraint plus six evolution equations for the six EM 
field components. One traditional approach to decouple the equations is to attempt to 
construct a second-order differential equation for each of Ei and Bi and it is shown in jB] 
that some of the 2+1 + 1 EM perturbations do satisfy decoupled equations. In this paper, 
an alternative method is described which proves to be successful in decoupling the entire 
system. We seek to decouple this system naturally by using the linear algebra techniques 
specified in |Appendix B| There exists a complex conjugate pair of combinations of Ei 
and Bi which will naturally decouple the system, and a new complex tensor is defined 
accordingly, 

q^.= E^ + iB^. (55) 
The system (J5T|) -(|51 jl is then given by, 

D k V k = ^e, (56) 

(C n - ~K)% + 2 Ai k V k + i e km {D k + h k )V m = - X - i { . (57) 

It is simple to show that the system for the complex conjugate (^j) can be found by 
taking the complex conjugate of (}56|) and (}57j) . 

The decoupled system for the EM fields can be determined entirely from either 
or Each system involves four complex equations governing three complex quantities. 
The complex nature of these quantities is arising due to the invariant structure of the 
source less equations under the simultaneous transformation (E^ — > B^B^ — > —E^). 
The use of a complex combination of the EM fields is well established throughout the 
literature and dates as far back as The purpose of this first stage of decoupling is 
to emphasise that the use of complex quantities is not a random choice, but instead a 
natural construction arising due to the inherent structure of the equations. 

The first stage decoupling has been completed in that now we have a system 
involving one tensorial quantity. The next challenge is to decouple the individual 
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components of The system for can be decomposed according to the 2 + 1 + 1 
decomposition, 

(C N -K)$ + d A <!>A=\e, (58) 
(C n - ~K + E)$ + i e AB d A <S> B = -ij, (59) 

Cn-^K-Ej^A-i e A B (C N + A)$ B + i e A B d B <£ = ~I Ai (60) 

where ^ = $ M + $-/V M , $ := iV a \l/ a , $ M := S^^a and the three-current has been 
decomposed as = TN^+I^. The parity reversed form of (JBOJ) is derived by contracting 
it with completely anti-symmetric the two-Levi-Civta tensor (which is defined such that 

(JC N + A)$ A - i e A B (c n - - Zj$ B - d A $ = i X -e A B T B . (61) 

Ultimately, it is desirable to find a decoupled equation for $ such that the wave operator 
is present. Therefore, it is required to know how the scalar wave operator appears in 
both 3+1 and 2+1+1 form[6j which are respectively, 

□0 : = V a V a <p 

= {-{C n ~K)C n +{D k + h k )D k )(t ) 

= (- (£„ - K) C n + (C N - K + A)C N + d a d a ) (f) (62) 

where is any scalar function. 

In the fully non-linear case, Maxwell's equations are heavily coupled to the ADM 
and GEM equations through extrinsic curvature and acceleration terms appearing in 
the equations, and through the energy-momentum tensor. We now discuss how they 
simplify under a first-order perturbation. 



3.2. EM Perturbations to Non-Vacuum LRS Class II Spacetime 

The EM perturbations to an LRS class II spacetime consists of both EM fields and 
charged sources, i.e. e, X and X A become first-order quantities. This paper 

assumes that these first-order EM fields and charged sources are gauge invariant*. 
Therefore, due to the Stewart- Walker lemma [2*T]. they vanish on the background. This 
consequently places a restriction on the background energy-momentum (p, P. J ', II) that 
it also be non-charged. 

It is also important to consider how these first-order EM fields and sources feedback 
to the geometry. For example, a typical energy-momentum tensor for these EM 
perturbations, which includes a simple interaction term, is |3J 

T^v = FfiaF a u — -g 9iu> FapF afi + g^vA a J a , (63) 

They are gauge invariant under infinitesimal coordinate transformations. 
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where A^ is a four-potential defined F^ v := 2V^A u y This clearly involves second-order 
terms at most and will therefore vanish when truncating at first-order. This simplifies 
the calculations as there are no first-order gravitational effects that need to be accounted 
for, and thus the geometry remains fixed. 



3.3. Generalized Regge- Wheeler Equation 

The outline for the derivation of the decoupled covariant equation for $ is given here. 
First take the Lie derivative with respect to n M of (JoTJj) and use PS]) to interchange 
the Lie derivative with respect to n M and two-derivative acting on $^4, then (160)1 is 
substituted for £ n $A- Then (Hfijl is used to interchange the Lie-derivative with respect 
to N^ 1 and the two-derivative when acting on $,4. Finally, (|HH)l and (jH^j) are substituted 
to eliminate any remaining $^ terms to give, 



(£ n + S - - K)£ n - (£ N + A - 2 K) £ N - d a d a -2K g + p- P -U + 2A 



= i X - e^d a I p - \ \{£ N + A - K) e + (C n - K) 1] . (64) 

This is the Regge- Wheeler (RW) equation[Tn] generalized to LRS class II spacetimes 
with a non-charged background fluid. There are actually two independent equations 
here, one for each of the real and imaginary components of $, which correspond to 
those developed in 6 J. 

By inspection of the differential operator above with ()62j) . there are clearly 
additional derivative terms left over once the true scalar wave operator has been 
identified. It is possible to rescale $ to obtain a true wave equation with a potential and 
a source defined in terms of energy-momentum quantities on the spacelike three-slices. 
This procedure is also presented in [Oj and the scaling is given by 

M := $ exp (Q) where C n Q := \y, - \ K and C N Q := --JC. (65) 

It is important to show that these equations for Q are in fact integrable, i.e. they must 
satisfy flU, 

(£ N + A)£ n n-(£ n -Z-]-K)£ N n = 0. (66) 

o 

which is achieved using the 2+1+1 field equations. Under this scaling, the usual scalar 
wave operator appears with additional potential terms arising, 

(□ - V)M = S, (67) 

jt The corresponding equations presented by 6 have slightly different numerical factors residing in the 
potential. The results presented here have been checked by expressing the equation using both the 
Schwarzschild metric (for which A, IC,£ ^ and T^,K — 0) and the Painleve Gullstrand metric (for 
which A = and K,, £, S, K 7^ 0) and the corresponding coordinate transformation yields consistent 
results. This cannot be achieved using |H]'s results. However, it is of primary importance that 6 's 
decoupling is not affected. 
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p-n + 2A, 



-i e aP d a l p + (£ N + A--JC)e+ (C n - K)X 



(68) 
(69) 



where 

s ■= ¥ 

In the absence of charged sources (S = 0), this wave equation may be separated out 
into a time/radial part and angular part. This is a because the potential is a function of 
t and r only. Such vacuum solutions can then be found using standard techniques; for 
example, spherical harmonic solutions arise for the angular components. These vacuum 
solutions offer invaluable information regarding how the EM waves interact with their 
background surroundings. However, in the presence of charged sources, this gives the 
exciting prospect of identifying the EM signature from a particular source. 



3.4- The Generalized Bardeen-Press Equations 

We now extend the results presented in [6] to show that the natural construction of 
complex quantities also yields a decoupled equation for $a- The derivation is analogous 
to how the RW equation was constructed. Take the Lie derivative with respect to 
of (jnOj) and use (|41j) to interchange the Lie derivative with respect to and the two- 
derivative when acting on $, and then substitute (|59|) to eliminate Also use ()44|) to 
interchange the two Lie derivatives when acting on $^ and subsequently, substitute (jfiUj) 
to eliminate £ n $A- Then use (I42|) to interchange the Lie derivative with respect to 
and the two-derivative when acting on $, and use (J58J) to eliminate £jv$. Finally, use 
(I60|) and its parity relationship (|61|) to eliminate any remaining two-derivatives acting 
on $, to show that 



[(C n -K)C n -{C N + A-K)C 



N 



d B d 



B 



+ 



K g -(£„ + £- -#)(£ 



$4" 

3 ; 



ie A B [(2A + K.)C n + 3XC N ]$ B 



(C N - K)A 



B 



(A 



K){A-\lC) 



{£n + A 



ire/ [d B l-(C 



N 



K)Xb] - 



(£ n + Z--K)l A -d A e 



(70) 



. 1 

This covariant complex equation has clearly decoupled from $, and the next difficulty 
is to show how to decouple the two independent components of 

This equation can be expanded in component form, by specifying an arbitrary LRS 
class II form of the metric, and then expressed in terms of large, untidy, matrices. The 
natural decoupling methodology used in Section I3~T1 and [Appendix B| can be employed 
again, and this again demands that complex combinations be formed to naturally 
decouple the equations. 

However, it is not necessary to specify the metric at all. Instead, a natural covariant 
decomposition of all two-surface quantities can be achieved by first defining a complex- 
conjugate pair of vectors (m /1 ,m /1 ) which satisfies the following relationships: 



m a m a 



m a m a 



0. 



m a m a 



0, 



r = 2m ( V ) . (71) 
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These complex-conjugate vectors are orthogonal to both n M and N^. Consequently, they 
can be raised and lowered using the two-metric. 

The remaining equations on the two-surface can now be irreducibly decomposed. 
Consider any two-tensors V^, V^ u (i.e. a contraction of any of their indices with n M or 
will vanish), where V^ v is symmetric and trace- free with respect to the two-surface 
(S a/3 V a p = 0). By using the properties specified in (J7TJ), they are both irreducibly 
decomposed according to 



V, 



{y a j3 , m a m^)m lJL m l , + (Vaf3m a m li )m f j i m 



(72) 
(73) 



and the generalization to tensors of different type is evident. Therefore, &a is irreducibly 
decomposed as 



$4 = M^mA + Mi 
where the scalars are defined M. ffi 



m A , 

= m A § A and 



(74) 



m A ^A- It is always possible to 



perform a rotation such that m M — > e~ lv5 m M , which, by using (fTTj) . clearly leaves the two- 
metric invariant. Under this rotation both .M© and M.® are not invariant. However, 
their reconstruction back into $,4 using (fTlj) is invariant. 

The two decoupled covariant equations arising naturally from the decomposition of 
(J7UJ) are, 

[£ n + 2 7 - K - q (2 A + K)}C n M® -[C N + 2X + A-JC + 3qi:] CnM® - d a d a M (B 
(£ n - K + - p(2A + !C) - (£ N + A - !C + \)\ - 3luZ + K g 



-2 X a d a M, 



-q{C n -K){A 



2 ; 



q (£n + A 



K)A 



- 3 K) +X 
- 1 S 



d a Xc 



(75) 



[C n + 2 7 - K + q {2A + K)\£ n M® - [C N + 2 A + A- K - 3 gE] C N M® - d a d a M t 



+2 x a d a M®+ [(£„-K + 7 )7 

+q (C n -K)(A 



p (2A + K) - (C N + A - /C + A) A - 3 to E + K g 
X -K) + q(C N + A- fc)(±£ + 2 -K)+ x + d a Xc 



-(£„ + E - |fcT)(E + ^K) - (C N - /C).A 



2 * 



(76) 



where the source has been decomposed as := S^rh^ + S^m^. The newly defined 
coefficients arising in these equations are related to various derivatives and combinations 
of m M and m M . They are given in |Appcndix A for LRS class II, and it should be noted 
that in most cases they are related to other geometric LRS class II scalars. 

The covariant decoupled equations, (fTKj) and (|TK|) . are the direct generalization of 
the Bardeen-Press equations [2 El EE3 UZj to non-vacuum LRS class II spacetimes. In 
the absence of charged sources (S 9 = S® = 0), the operators are separable which is due 
to the fact that the LRS class II scalars are functions of t and r only and the two-metric 
is expressed in terms of the complex-conjugate vectors. Such vacuum solutions may be 
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obtained using standard techniques. In this case though, the angular solutions will be 
in terms of spin-weighted spherical harmonicsf^j. Furthermore, it is possible to specify 
energy-momentum on initial three-slices as the sources to the EM radiation. 

3.5. Example: Static Schwarzschild Background 

We now discuss the static Schwarzschild background to show that (ff5|) and (ff6|) are a 
generalisation of the BP equations. The static Schwarzschild line element is 

ds 2 = - (l - dt 2 + + r 2 (d9 2 + sin 2 Qd<\> 2 ) , (77) 

and the 2+1+1 background equations reduce to 

(C N + A-)C)A = 0, (C n - A - -K)K = and £ = AlC. (78) 
Therefore, the only non-zero components of the LRS class II scalars are, 



„ M ( 2M\* 2/2M , „ 2M 
A=-? 1 , K = — a/1 and £ = =-. 79 

Furthermore, the complex-conjugate vectors defining the two-metric is defined 

1 



_ (O,O,l,icosec0), (80) 
y2r 

but this is not a unique choice. Before proceeding, it is noted that the decoupled 
equations were checked by comparing them with those derived using the Newman- 
Penrose formalism [8]. The null vectors used for this purpose were defined according 
to [S] and they are also not a unique choice. Therefore, to put the equations in the 
exact same form as those presented in [Hj, one final scaling is requiredff , 

Me:=/#ffi and M®:=jB®, (81) 

where C^^f = A and C n f = d^f = 0. Finally, including the RW equation, we have 
three decoupled equations, 

[C n C n -{C N + A-2K)C N - d a d a - V (RW) ] $ = 0, (82) 
[{C n -2A-K)C n - (C N + 3A-2)C)C N -d a d a -2x a d a -V (B \B (B = 0, (83) 
[(Cn + 2A + K)C n - (C N - A-2K)C N - d a d a + 2 x a d a - Vy = 0. (84) 



Here the Lie derivatives reduce to C n = and Cn = V -L 11 ^:, and the potentials are 
V {KW) :=2K g , (85) 
V 9 :=~(C N + A-K)K:-x = K g - X , (86) 

V 9 := ~\(Cn -3A-lC)lC-x = K g + 2S-x. (87) 

ff In general the scaling is not necessary. The equations were checked by expanding everything in terms 
of coordinates using Maple 9.5 and they correspond exactly to those derived using the Newman-Penrose 
formalism [HJ. 
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The RW is given by (J82)) for the quantity $ which has zero spin- weight. Furthermore, 
(JHHj) and are precisely the BP equations El EH \U\ for the spin- weighted quantities 
of +1 and —1 respectively. 

4. Discussion 

We have shown how to derive a fully decoupled set of three covariant equations, 
which completely determines the EM perturbations about a non-vacuum LRS class II 
spacetime. One of these is the generalized Regge- Wheeler equation, whereas the other 
two are the generalization of the Bardeen-Press equations. The setting uses the 2+1 + 1 
formalism, and consequently, the perturbations are rich with physical significance. The 
equations can be expressed in terms of quantities such as the Gaussian and extrinsic 
curvature of two-dimensional surfaces thereby rendering the interpretation relatively 
simple. In vacuum perturbations, all the differential operators are separable and 
solutions can be found using separation of variables, and the usual spherical harmonics 
and spin-weighted spherical harmonics naturally arise. Furthermore, the sources to 
the equations are written in terms of charges and currents specified on initial spacelike 
three-slices. This gives the exciting prospect of modelling EM radiation from a particular 
astrophysical source, rather than the more abstract vacuum pertubations. 
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Appendix A. Complex-Conjugate Identities 

In LRS class II the non-trivial identities and definitions are: 

7 := m a C n rh a , 7 + 7 = £ - ^K, C n m a C n rh a = -j 2 , 

\:=m a £]yma, A + A = — /C, H^m a L^m a = — A 2 , 

Xa ■■= m a d^m a , x-={d a m l3 )(d a m () ), p := i m a e aP C n m p , 

q : = im a %e a/3 = ±1 and u> := im a e al3 CNfnp- 

In the static Schwarzschild spacetime, with the complex-conjugate vector defined by 
flHUj) . they reduce to 



(A.l) 



7 = 7 = p = p = d a Xa = 0, A = A = u = to = -^/C, 

cot 2 # 

q = l, Xy. = [0,0,0,1 cos 6\ and x = — — ■ 



(A.2) 





Appendix B. Linear Algebra: Decoupling 

Consider the system given by 

L 1 E + L 2 B = ® and L X B - L 2 E = 0, (B.l) 

where L\ and L 2 represent differential operators and E and B are any scalar fields. This 
system has the property that it is invariant under the simultaneous transformation of 
E — > B and B — » — E. This system can be expressed in a matrix form as 

(B.2) 

where M :— ( ^ J J is the matrix responsible for coupling E to B. M can be 

written in terms of its eigenvalues, diag(D), and corresponding eigenvectors, col(P), 

i -i \ ( -i \ x ( -i 1 

1 1 y \ i y 2 v i 1 

is diagonal, it is clear that by multiplying (jB.2|) by —2 i results in the decoupled 

system, 

Lx{E + iB) +i L 2 (E + iB) = and L^E — \B) — i L 2 (-E — i 5) = 0. (B.3) 

Thus a complex-conjugate pair of equations arise. This result can be generalised to 
tensors of any type without loss of generality provided the invariance is satisfied. The 
matrix M needs to be written in block form with blocks of zeros or the identity matrix 
to compenstate for the number of dimensions. 



according to M = P D P' 1 = ( i )( n .)|( . 1 ). Therefore, since D 



